Estimating multipartite entanglement measures 



Andreas Osterloh 1 ' 2 and Philipp Hyllus 1,3 

1 Institut fur theoretische Physik, Leibniz Universitat Hannover, Appelstr. 2, 30167 Hannover, Germany. 
2 Fakultat fiir Physik, Campus Duisburg, Universitat Duisburg-Essen, Lotharstr. 1, 47048 Duisburg, Germany. 
3 BEC-INFM, Dipartimento di Fisica, Universita di Trento, Via Sommarive 14, 1-38050 Povo, Italy. 

We investigate the lower bound obtained from experimental data of a quantum state p, as proposed indepen- 
dently by Giihne et al. and Eisert et al. and apply it to mixed states of three qubits. The measure we consider 
is the convex-roof extended three-tangle. Our findings highlight an intimate relation to lower bounds obtained 
recently from so-called characteristic curves of a given entanglement measure. We apply the bounds to esti- 
mate the three-tangle present in recently performed experiments aimed at producing a three-qubit GHZ state. A 
non- vanishing lower bound is obtained if the GHZ-fidelity of the produced states is larger than 3/4. 
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I. INTRODUCTION 



Since entanglement has been recognized as a possible valu- 
able resource for quantum information processing fin, its 
analysis, detection and quantification are major goals IJ-Ql]. 
Except for simple measures that are polynomial invariants of 
degree 2, as e.g. the concurrence, already the calculation 
of that measure for a known mixed state poses a hard prob- 
lem. In the laboratory, however, the states to work with are 
mixed states with major weight on a desired entangled state 
and additional uncontrolled admixture of different states due 
to systematic or nonsystematic errors in the preparation pro- 
cess. Even though in priciple full state tomography can be 
performed, it is experimentally increasingly expensive with 
the number of qubits. 

It is for this reason that entanglement witnesses J§| play 
an important role for experimental detection of entanglement. 
Witness operators are constructed such to have negative ex- 
pectation values only on states that carry a specific class of 
entanglement; all states that do not belong to this class have 
positive expectation value of the witness. Then, a negative 
expectation value implies that the mixed state p carries the 
specific entanglement class detected by the witness. Recently, 
also quantitative estimates have been obtained for a variety 

of entanglement measures from (ejj. experimental) values 

for one or more witness operators |7 



10 



, and for non-linear 
1311 . by using meth- 



combinations of expectation values 1JJJ 
ods of convex optimization lfl4ll . Further lower bounds have 
been obtained using a different approaches in a similar con- 
text lfl5l - l20ll . A connection of the problem of entanglement 
estimation from uncomplete information to Jaynes principle 
has been discussed already in |2lll . 

Here, we investigate the tight lower bound for convex-roof 
extended measures as introduced in Refs. |0 Htl for mea- 
sures quantifying true multipartite entanglement in order to 
compare them with lower bounds obtained from different 
premises l22ll . The computation of the lower bound involves 
a supremum over some parameters and an infimum over pure 
states. One of our central observations is that the bound can be 
alternatively computed in a way which involves an infimum 
over pure states only, followed by the convexification of the 
resulting function. This observation establishes a connection 



between the methods from Refs. and Ref. iB . 

To this end we focus on the three-tangle, where a specific 
analytic solution of the convex roof can be used as a bench- 
mark. 

The work is laid out as follows. In the next section we 
briefly sketch the estimation method from witness operators 
as proposed in Ref. |0 [Ml an d introduce the three-tangle and 
general multipartite entanglement witnesses. We then turn to 
the application of the method in Section [Hi] using expecta- 
tion values of entanglement witnesses designed to detect true 
three-partite entanglement. We first consider a specific class 
of rank two mixed three qubit states for two different witness 
operators. The estimation using the data of the first witness is 
directly related to the results of Ref. i22tl . Then, we compute 
a lower bound for an important witness, which can directly 
be used to estimate the three-tangle produced in recent exper- 
iments. Finally, we report on further interesting observations 
as the effect of common symmetries of the Witness operators 
and the entanglement measure in Section IIV1 where we also 
sketch a proof for the possible reduction to the related problem 
for pure states. Section [Vl contains our conclusions. 



II. BASIC CONCEPTS 

A. Lower bound on entanglement 

We consider the following situation: assume that a source 
can be described by the (unknown) density matrix p and that 
K expectation values Wk = tr [pWfc] are measured and then 
collected in a vector w. The operators Wk are further assumed 
to be witness operators. The lowest value of the entanglement 
of the state p consistent with the measurement results is for- 
mally given by the solution of the problem 



inf E{p) 

p 



tr [W k p]=w k 



(1) 



In Refs J3,0] it has been shown that the solution of this opti- 
mization problem can be approximated from below with affine 

w - E(J2k r kWk)), where 



functions by e(w) = sup r 
E(W) = sup p (tr [pW] — E(p)) is the Legendre transform 



2 



of E. It was further shown that for pure state entanglement 
measures extended via the convex roof construction as 

E(p):= irf (2) 
p=Y, i Pi\ipi)(ipi\'~ 

where £\ p t = 1 and p. t > 0, the Legendre transform can be 
computed by optimizing over pure states only, leading to 

K 

e(w) = supinf ( V r k (w k ~ W k |V>» + £(|V») . (3) 
r \4>) \ k=1 

This bound is tight due to the convexity of the problem 0, 

A key observation is that Eq. (f3j) is the dual problem to the 
minimization of E(ip) on pure states subject to the witness 
conditions. It is therefore solved by an approximation of the 
pure state problem 



inf E(ip)\ 



(4) 



from below with affine functions Hence the result coincides 
with the function convex hull of the pure state problem l37ll . 
In other words, in order to solve the problem ([T]) for convex- 
roof extended entanglement measures, we could alternatively 
solve the problem (01 and convexify the resulting function. 
In cases where K is small, this could reduce the computa- 
tional cost. This point is illustrated in Section IIII Al where 
we explicitly solve (01 and (01 for a restricted situation, and 
the equivalence of both problems is highlighted explicitly. 
The origin of this equivalence is further investigated in Sec- 
tion llV Al for arbitrary K. 

The entanglement measure we consider is the three-tangle 
and its square; it distinguishes the two classes of global entan- 
glement for three qubits: W and GHZ. 



B. Three-tangle 

The three-tangle of a general pure three qubit state ex- 
panded in a product basis = k=o Ni^) * s gi ven 
by Q 

t 3 = 4|di-2da + 4d3| (5) 

d i = V'oooV'iii + V>oWiio + ^(Wioi + ^iooV&i 
d-2 = 0oooV j iiiV'oiiV ; ioo + 0oooV , mV'ioi'0oio 
+^000^111^110^001 + V'onV'iooV'ioiV'oio 
+^011^100^110^001 + ^ioiV'oioV'iioV'ooi 
^3 = V'oooV'iioV'ioiV'oii + V'mV'ooiV'oioV'ioo- 

It can be extended to mixed states via the convex-roof con- 
struction (0). 

The three-tangle is non-vanishing on true 3-partite entan- 
gled states only and vanishes on any bi-separable state such as 
|^ab)<8>|0c) for the parties A, B, and C, but also forW states 
introduced below. Also the convex-roof extended measure 
has this property for mixtures of bi-separable states, possibily 
for different partitions of the parties. Summarizing, it dis- 
tinguishes the two classes of three-partite entanglement that 



are inequivalent under stochastic local operations and classi- 
cal communications (SLOCC) 1 24] . The representative of one 
class is the W-state 



|W) 



1 



V3 



(|001) + |010) + |100)) 



(6) 



and the representative of the other class is the GHZ state |25| 



|GHZ> = i=(|000) + |lll>). 



(7) 



The three-tangle vanishes on all states of the W-class while 
it is non-vanishing for states of the GHZ-class. This classi- 
fication can be extrapolated to mixed states in the sense that 
a mixed state belongs to the W-class if it has a decomposi- 
tion into pure states exclusively out of the W-class 12611 . In 
analogy, a mixed state belongs to the GHZ-class if it has a 
decomposition into pure states of the GHZ-class. 



C. Multipartite entanglement witnesses 

Multipartite entanglement witnesses can be constructed 
as JH 



W = afl- 10X^1, 



(8) 



where a = max| Zero ) |(0|Zero)| 2 is the maximal overlap of 
|0) with any state for which a given entanglement measure 
£7(|Zero)) = 0. In the case that global entanglement is con- 
cerned, |Zero) would be all biseparable states, whereas also 
VK-states were included if only states with a GHZ component 
are of interest IH [H, |H] . We will highlight the latter situ- 
ation and use the three-tangle, E := t$, since it is non-zero 
only for states of the GHZ class. 

For such witnesses, Eq. (O becomes 



K 



K 



sup I - J> fe (0,1 p |0 fe ) + inf (£(|0)) + ]Tr fe |(0 fc |0)| 2 ) 

r I fc=l k=l 

(9) 

This tells us that the trivial part ai of the witness is irrelevant 
for the estimation. Hence, it does not seem to be essential that 
the expectation values are measured with respect to witnesses; 
an operator which projects onto a suitable subspaces seems to 
be sufficient J^]. 

For a single witness of this form it is straight forward to 
check that r > leads to e(w) = 0; indeed, in this case the 
infimum over all states |-0) would be those states perpendicu- 
lar to |0) for which E(ip) = 0. This is always possible for a 
single witness. Another peculiar case is when (0| p |0) = 0. 
In this case the estimate is independent of the density matrix 
under consideration and hence can only be zero. This cor- 
responds to a bad choice for the witness operator. For the 
three-tangle, choosing |0) outside the GHZ class also leads to 
e(w) = since then the infimum in (O is reached by choosing 
|V') = |0) forr < 0. 

The above reasoning has interesting implications. First, r 
can be taken to be negative for a single witness operator de- 
tecting the class of entanglement measured by E. Second, it 



will admit to relax the restriction for W to being a witness, 
and a wider class of observables might lead to a reasonable 
estimation of E. In these cases, r can possibly also be pos- 
itive i38ll . For more than one witness, the situation is more 
complicated and then even observables that are not useful at 
all for detecting entanglement could possibly be useful for im- 
proving the lower bound. 



III. APPLICATION 

We illustrate the results of the previous Section by explic- 
itly solving the optimization problem <f3j for special cases 
where knowledge about the state is assumed to be given and 
by showing explicitly that the result coincides with the one ob- 
tained by solving problem and convexifying the solution. 
Then, we numerically solve the general problem for a general 
witness which allows to estimate the three-tangle from exper- 
iments where the GHZ state was produced and the fidelity 
(GHZ | pcxp | GHZ) was estimated for the produced state p exp - 




FIG. 1: (Color online) The graph shows how T$(q, 0) varies with q. 
The graph reaches at q = and go ~ 0.627. The red dotted curve 
is the convex-roof T3 (p) for the mixed state d 10b - Inset: Deviation 
of T3(q, 0) from the straight line going from (qo, T3(qo)) to (1, 1) in 
the interval [qo, 1]. Note that the behaviour of Ts(q, 0) changes from 
convex to concave. 



A. Restricted optimization 

We will demonstrate the Legendre transform at work on 
a simple example, where the Hilbert space for minimization 
with respect to ip is restricted to the range of the actual den- 
sity matrix. This amounts to an idealized situation where the 
experimenter knows precisely which states are possibly pro- 
duced in the laboratory to be part of the resulting mixed state. 
The specific example is about the three-tangle of the three 
qubit mixed states 



p(p) = p 7TGHZ + (1 - P) 7TW , 



(10) 



for which the analytic convex roof is known 12811 . Here, 7Tghz, 
7Tw are the projectors onto the |GHZ) and |W) state, respec- 
tively. The vectors in the range of p(p) are given by 



\Z(qA)) = V^IGHZ) 
with its three-tangle given by 



J4> 



s/T^q\W) (11) 



8V6 
~9~ 



3i<j> 



T~z{qA) := T 3 (Z(q,(f>)) = 

112) 

In the following, we want to apply the method of Refs [[Zllll] in 
order to get lower bounds for the three-tangle of general states 
g in the GHZ-W subspace, using two qualitatively different 
single witnesses by performing the optimization of Eq. (0). 
We the compare the result with that obtained by solving the 
problem (0]i and by convexifying the solution. 



1. GHZ witness 
The first witness we would like to consider is 

WgHZ = at — 7TGHZ- 



For a = this is a witness for multipartite entanglement, 
while for a = |, it detects states of the GHZ class only l26ll . 
For this witness, w = tr [Wg] = \ — p, where we defined 
P = (01 Q 10)- We perform the infimum with respect to \i[>} 
only over the restricted set of states \Z(q, <j))) in the range of 
g. This is useful for obtaining a bound on the tangle of the 
rank 2 state g and relies on the assumption that we know the 
subspace that g lives in. As we will show below, knowing the 
expectation value of the witness dT3b corresponds to knowing 
q in Eq. ( fTTT >. Hence the only free parameter in the pure state 
problem © is the phase <f>, which is minimized for = 
1281 . By the argument of Section [n] the lower bound to the 
problem is then given by the function convex hull Ts(q) of 
T3(g,0) (see Eq. (fT2b and Fig.Q]), corresponding to the con- 
vexified solution of problem Q. It has emerged as so-called 
characteristic curve from a different approach to obtain lower 
bounds on entanglement measures pursued in Ref. M22I1 . 

Let us compare this result to the problem problem ® which 
can be solved explicitly in this simple case. We have to deter- 



e(w) = supinf 

r q,<f> 



(r(q-p)+T 3 (q,<l>)\ 



(14) 



(13) 



and we show in what follows that e(w) = fj,(p) is obtained 
also in this case. As mentioned above, the infimum over <j> is 
always obtained for <f> = so that we have to optimize over 
q only. T^{q, 0) is plotted in Fig.[T] Graphically, optimization 
in Eq. ( fT4l means to tilt the x-axis about the fixpoint (x, y) = 
(p, 0), resulting in a curve with local minima for each tilting 
slope r. The largest of these minima is the estimated lower 
bound for the entanglement measure, here T3. This seesaw 
argument is shown for two cases in Fig. [3] 

As already discussed earlier, assuming r > leads to an in- 
fimum at q = for allp £ [0, 1] and we are left with the trivial 
bound e(w) = sup r (— rp) = 0. So we take r < and focus 
\r\q + T3(<7, 0)), formally looking for 



at the infimum inf „ 



0.6- 



0.4- 



0.2- 



< 



0.5 





FIG. 2: (Color online) The estimate for from a diagonal wit- 
ness (red curve) together with the square of the corresponding char- 
acteristic curve (black). It can be seen only in the inset that very 
close to p — 1 this characteristic curve still fails to be convex; the 
lower bound correctly gives the function convex hull. The small inset 
zooms into the region, where the characteristic curve fails to be con- 
vex. It shows the numerical error, which is two orders of magnitude 
below this deviation. 



local extrema in intervals where r 3 (q) is twice differentiable. 
The differential minimum condition is then 



M = Ts(g,0) 
h(q,0) > 0. 



(15) 
(16) 



The second condition simply demands that T 3 (g,0) be con- 
vex (see the discussion in Section ITVAI for the general case). 
Where T%{q, 0) is concave, the infimum is found at the bor- 
der of some interval, which has to be determined. We have 
to find the global infimum and three cases have to be treated 
separately. 

Unregarded the value of p there are three separate regions 
for the variable r: (i) > r > ro, (ii) ro > r > n, and 
(iii) j"i > r with ro, ri to be specified below. The three cases 
can be understood graphically from Fig. Q] In region (i), the 
infimum is taken at qo = 3 ^J|^ ~ 0.627 until the slope 
of the function —\r\q + r 3 (<?, 0) (for q > qo) at q = qo van- 
ishes. This happens for ro = — r 3 (go,0) ~ —2.52. In region 
(ii), the function — \r\q + r 3 (q, 0) has a unique local minimum 
at some position q r € [qo, q{\ given by Eq. ( fT5l ). where q\ is 
determined below. This local minimum is also the global min- 
imum. We see that — [r[q + Ts(q, 0) is convex in this region, 
when zooming into the plot, as shown in the right panel of 
Fig-Q] Finally, n is the value for r below which the infimum 
is located at q = 1. The values of r\ and qi can be computed 
as follows: at r = r\ the infimum is reached both at qi and at 
9 = 1. We therefore have r\q\ + T 3 (qi,0) = r% + r 3 (l,0), 



and hence r\ 



l— T&(gl) 



Also at this point, the slope of 



1-91 

r 3 (q, 0) equals that of the straight line connecting the points 
(<Zij T 3(9i; 0)) an d (1) I)- hence r\ = —73(91). Both condi- 
tions determine qi = i + ^1^-^/465 ~ 0.70868. Reinserting 
p, we are left with the following optimization problems in the 
three regions: (i) sup r r(g —p), (ii) sup r r{q r —p) +T 3 (q r , 0) 



0.7 -„ 



OO 



0.6 



3 

+ 

S 0.5h'' 
e 0.4 



r=-2. 19799 
r=-2. 69799 
r=-3.19799 



0.7 



0.75 



0.8 



0.85 



0.9 0.95 



FIG. 3: (Color online) Here, we show the seesaw argument for the 
three-tangle. The upper panel shows the situation for p = 0.5 (cor- 
responding to region (i)). Besides the optimal value r = (green 
full line), one positive and one negative r value has been considerd; 
the largest absolute minimum of the tilted function is clearly zero. 
The lower panel considers p — 0.85 (corresponding to region (iii)). 
Here the optimal value r = n is shown (green full curve) next to 
two close but deviating values. 



(qr € [<7oj<7i]X while in region (iii), the infimum is given 
by — |r| + 1, and e(w) = sup r (— |r|(l — p) + 1) = 1 — 
I ri I ( 1 — p) . This curve is the straight line connecting the points 
{qii t 3(qi)) an( l (1) 1)' as already mentioned before. 

The question we address now is, which of the above de- 
picted three regions corresponds to a given p. For p < qo, 
choosing r from region (i) yields e(p) = since r < and 
qo — p > 0. In regions (ii) and (iii) the infimum is at negative 
values (see Fig. [3}. For qo < p < qi, we find e(j) = ro(p — qo) 
in region (i), and q r = p in region (ii), hence = r 3 (p, 0). 

In region (iii) the infimum is at negative values again. Due 
to the convexity of r 3 (p, 0) in p g [q , qi], r 3 (p, 0) > r (p - 
qo), and hence e(w) = r 3 (j5, 0). Finally, for qi < p < 1 the 
supremum of the infima is located in region (iii), and e(w) = 
1— |ri|(l— p). Altogether, this yields e(w) = f 3 (g) as claimed 
above. 

This is further illustrated by numerically applying both 
methods to r| ll39ll . as shown in Fig. [2j The region, where the 



5 



0.5 



We obtain 




FIG. 4: (Color online) For visualizing how the estimated lower 
bound is decreased by the off-diagonal part of the witness J 1 7b . we 
assume the term yj q(l — q) cos <f> to appear without r as a prefac- 
tor for a moment. It then simply modifies the characteristic curve as 
shown here for uj = 1. The function convex hull of the lowest en- 
velopping curve (thick full red line) would then constitute the lower 
bound estimate from the off-diagonal witness. However, this esti- 
mate would still lead to the correct value 1 at p = 1. Even this is 
spoiled when the contribution from the off-diagonal part of the wit- 
ness is actually r-dependent, as is the case here. 



function convex hull must be applied is clearly distinguished 
from the numerical uncertainty. 

Quite generally, lower bounds can be obtained from the 
characteristic curve of a given entanglement measure l22ll . 
Here this lower bound coincides with that bound obtained 
from the approach in Refs and even gives the analytic 

convex roof. In this sense, the bound is tight. 



e(p) = supinf (r(q — p + 2^q(l 

r 9,0 V 



q)Re uje^) + T 3 (q,, 



(18) 

In this case, the infimum cannot be taken separatly for q and (f>. 
It is clear from the specific situation that the off-diagonal part 
would result in a decrease of the lower estimate of the three- 
tangle (see Fig. |4j». The fact that r is a prefactor of -J q(l — q) 
in the function to be optimized leads to the feature that even 
the estimate at p = 1 is below the exact value 1. We an- 
alyze this case in detail because it exhibits relevant features 
an improper choice for the witness might have. This type of 
experimental "error" might be due to e.g. incomplete knowl- 
edge about relative phases in the state under consideration or 
a not perfectly symmetric setup for the production of e.g. a 
GHZ state. The restriction to a two-dimensional Hilbert space 
admits to depict both possible effects and their origin. 

Also in this case we determine the bound e(w) with the two 
methods. We first solve problem (|3}, consisting in numerically 
reaching the supremum of the infimum. Then, we apply the 
witness restrictions directly to the pure states (see Eq. (01)). 
Here, this leads to 



p + 2u> 2 cos 2 4> ± 2|w cos (j>\\/ijj 2 cos 2 <fi + P P 2 

qmin{P,(p) = : , , 2 J"- 

1 + ito z cos z <p 

(19) 

where cj cos <f> < must be imposed, and the infimum is then 
obtained for the minus sign in ( fT9l . 

In the second case, the lower bound of the three- 
tangle is then obtained as the function convex hull of 
miri0 T3(q m i n (p, <fi),4>)- For r|, this curve is indicated with 
circles in the left panel of Fig. [5] together with the result from 
the numerical approach; their agreement is perfect [39]. 



2. Off-diagonal witness 

Clearly, the previous discussion marks an ideal situation 
in that the infimum over the admissible states 4> was easy 
to handle. As already mentioned in Section [Til it is essen- 
tial for the operator W to have a nontrivial overlap with the 
set of density matrices of interest. With nontrivial we mean 
that tr p W 7^ const. For GHZ-W mixtures it implies that 
W must be able to distinguish both states. Consequently, 
W = -a \W) (W\ - /3 |GHZ) (GHZ| can be employed (un- 
less a = (3), but also W" = - | 111) (111 |, both without 
changing the resulting lower bound as obtained above for the 
particular GHZ witness Wghz = 5! — ttghz- We emphasize 
however that this is due to the restricted Hilbert space we con- 
sider here. To point it differently: more ab initio information 
about the state admits more freedom to the observables in or- 
der to get sensible lower bounds. 

A more generic situation occurs when we admit for off- 
diagonal operators like 

W skcw = -^ghz - u |GHZ) (W\ - \W) <GHZ| . (17) 



B. Unrestricted optimization 

Now let us consider the situation relevant for most experi- 
ments. A state q is created and w = tr [pWghz] is measured 
ll40ll . From this value we would like to obtain an estimate of 
the tangle. This corresponds to obtaining a bound on the three- 
tangle given the fidelity of the state with respect to the GHZ- 
state, p = (GHZ| p |GHZ), since w = a — p. We have solved 
the problem (0) for this case numerically. The results are plot- 
ted in Fig. [6] In Table I, the corresponding bound on the three- 
tangle are listed for experiments where p has been measured. 
It is a curious fact that the results agree with those obtained 
by solving d4]i, where the optimization was performed only 
over states which are symmetric under the exchange of par- 
ticles. We add some more comments on this curious fact in 
Section HVBl 

The problem restricted on pure states amounts to finding 
the steepest descent of the three-tangle when the GHZ state is 
superposed with some other state out of the orthogonal com- 
plement, which is numerically found to be reached using only 
the states W and its bit-flipped form, W. The value zero is 
assumed for the weight p = 3/4 of the GHZ state. For GHZ- 
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FIG. 5: (Color online) The estimate for r| from an off-diagonal wit- 
ness for various relative weights ui := a/7 of the off-diagonal part. 
The insets show the deviation Ar| to the exact convex roof for rf. 
Increasing uj the estimate quickly decreases and finally becomes the 
trivial bound, which is reached well before u = 00, which means a 
purely off-diagonal witness. This is only due to the special choice of 
basis we did: the off-diagonal part does not detect neither GHZ nor 
W state. The lower panel shows that negative uj is more destructive. 
The circles and squares indicate the solution obtained from setting 
the part proportional to r of the functional to zero; they lie on top of 
the numerical curves. 

states mixed with white noise 

77TGHZ + (1-7)1/8, [0,1], (20) 

which we will refer to as noisy GHZ states, this leads to zero 
three-tangle at 7 = 5/7. This coincides of course with the 
known bound for GHZ-entanglement in these states: by con- 
struction, if p > 3/4 they are detected by the optimal GHZ 
witness Wq^ z = |l — 7tghz ( see Ref- It26l0 . Since this is 
a lower bound, and we do not know whether this bound is 
tight for the noisy GHZ states, we give an upper bound for its 
three-tangle atp = 3/4. 

To this end, a decomposition of the noisy GHZ state must 
be constructed, and we use our result that the minimum value 
for the three-tangle is assumed for symmetric states contain- 



FIG. 6: (Color online) The unrestricted optimized lower bound for 
the three-tangle tz is shown for a given GHZ fidelity p. This captures 
also the GHZ-identity mixture. Two different methods are compared: 
full scale numerical optimization and the optimization restricted to 
symmetric states satisfying the constraint given by the GHZ wit- 
ness 1 11 — 7tghz, which fixes the value of p. Both methods lead 
to the same result (modulo taking the function convex hull in the lat- 
ter method), as it should be. Inset: close to p — 1 the curve obtained 
by implying the witness constraint directly on the pure states is seen 
to lie above the curve resulting from numerical optimization; in this 
region, the pure state curve is not convex. The function convex hull 
agrees with the curve from the full scale optimization. 



Experiment 


GHZ-fidelity p 


<P) 


ion trap f291 


0.86 ±0.03 


0.42 ±0.12 


diamond centers [301 


0.87 ±0.06 


0.46 ±0.24 


ion trap [311 


0.979 ± 0.002 


0.914 ± 0.009 



TABLE I: (Color online) Lower bounds on T3 for recent experiments 
where the estimated p exceeded 3/4. The used data is the same as 
that used for Fig. [6] 



ing W and W alone and with equal weight. Imprinting local 
phases, we can induce relative phases for the W and the W 
state, that are mutually inverse. When keeping the phase of 
the GHZ state fixed and averaging over these relative phases, 
this leads almost to a state of the form d2-0T >. but with the or- 
thogonal GHZ state with relative minus sign, GHZ_, missing. 
Convexity of the convex-roof applied to the resulting decom- 
position leads to an upper bound of 1/9 for the three-tangle of 
the noisy GHZ state (|2Q| with GHZ weight p = 3/4. 

An interesting question to ask here is for the full convex set 
of states with zero three-tangle and GHZ weight 3/4. If the 
interior of this polytope contained a state of the form d20l ), 
this would mean that the lower bound was even tight for 
these states. Otherwise the bound is not tight and the thresh- 
old weight for GHZ-entanglement (as measured by the three- 
tangle) would definitely be smaller than 3/4. 

A lower bound p su b = 1/3 and 7 su b = 5/21 for the weight 
of GHZ in a state (l20b without three-tangle can be obtained 



7 



from a phase average of states of the form 

1 - 7 

7(pottghz + (1 -PoKwJ H — (ttghz_ +7T Wq ) 

which are known to have zero three-tangle [32]. The index 
a indicates that relative phases have been introduced into the 
state by means of local phases. 

C. The effect of more witness constraints 

If the expectation values of more than one witness opera- 
tor has been measured, then the class of states over which the 
optimization is performed is restricted. For instance, the prob- 
lem discussed in Section lTlI Al is recovered if p = (ttghz) and 
P = (t£w) are known and in addition p + p = 1 holds. In 
Ref. 12711 . p and p have been measured with a setup intended 
to produce a W-state. As expected, the lower bound on T3 
obtained from the data is equal to zero. 

Note that if the pure state problem © returns zero as an 
estimate, the final lower bound is also zero, since the function 
convex hull is smaller or equal to the original function. 

IV. FURTHER OBSERVATIONS 
A. Restriction to pure states 
The problem <£3j can be written as 

supinf {r(w- W{x)) + E{x)} , (21) 

where x is assumed to be a minimal set of parameters that 
uniquely describes the pure states of the system. It is worth 
mentioning that due to VfV?{. . . } =0, the order of the ex- 
trema is relevant. Therefore, we consider at first Vg{. . . } = 

0, which leads to r ■ W'{x) — E'{x). If W' is invertible in a 
neighborhood of the optimal point this is formally equivalent 
to 

r = V #J B, (22) 

that is, the vector of Lagrange multipliers equals the gradient 
of the entanglement measure with respect to the witness val- 
ues. The consecutive supremum in r, inserting this condition, 
leads to 

W(x) = w, (23) 

1. e., the pure states themselves already satisfy the experimen- 
tal witness contraint. 

A remark is in order: the number of witness constraints will 
typically be significantly smaller than the number of parame- 
ters describing the states. Then, for each set w of witness 
constraints, there will be a corresponding submanifold in the 
parameter space x. In order that W' can be invertible, we 
have to take the infimum of E in that manifold and call it 
E c (w), which is the characteristic value of E as introduced 



in l22ll . and we can henceforth use w as the parameters de- 
scribing those states with minimal value for the entanglement 
measure subject to the contraint (ip\\V\ip) = w. 

For positive semidefinite second derivative we have an infi- 
mum; after straight forward algebra, this leads to the condition 
that the matrix with entries 

d 2 E c dE c 8 2 w k = d 2 E c 
dwidw m dw k dwidwm dwidwm 

be positive definite. This requires convexity of E c . The equal- 
ity holds where W' is invertible and the witness constraints are 
independent, i.e. we assumed d Wl u>k = Ski- Wherever E c (w) 
is not convex, the infimum will be assumed at the boundary 
of some interval, and the function convex hull will have to be 
taken at the end. 

Interestingly, invertibility of W' and g^-^g , when in- 
serted into the second derivative w.r.t. r leads to 

Vf>V P (. ■ ■ )l ex trema conditions = ~ (V% E c(&j) , 

which should be negative semidefinite, since we maximize 
with respect to r. This is consistent with our above result that 
V^Edw) had to be positive definite. 

This analysis generalizes the explicit calculation for a sin- 
gle witness performed in Section lTH A II It underpins the deep 
connection to the concept of characteristic curves as proposed 
in Ref. [22]. In the presence of non-convex regions in E c (w) 
the extrema are assumed at the boundary dl^ of some interval 
lis in u>-space; in complete an alogy to the discussions on op- 
timal decompositions in Refs l22l l28l l32ll this means that no 
pure state satisfying the witness conditions leads to the lower 
bound for this case. The lower bound is then achieved by 
mixed states made of some pure states in dl$, and the lower 
bound is affine within the whole interval 1$. For a single wit- 
ness, this boundary consists of precisely two states and then 
prescribes a mixed state solution of the problem. 



B. Simplification due to symmetry 

The problem also simplifies, when both the entanglement 
measure to be estimated and the observables measured in the 
experiment have a common symmetry [41]. Let us assume 
that both E(p) and tr pWk are invariant under p — > Qj(p) = 
QjpQ], where the operators Qj form a group Q with J ele- 
ments. The symmetry group encountered later on is the sym- 
metric group S n of qubit permutations; therefore we assume 
that J is finite in what follows. The result is more general 
and also applies to compact continuous groups. Let us define 
p = j J2j=i Qj (p) gi yen a state p. Clearly, Qj (p) = p holds, 
and then convexity and invariance of E with respect to J im- 
ply E(p) < E(p). Note that also p fulfills the experimental 
constraints Wk = tr WkP since also the witnesses Wk are as- 
sumed to be Q invariant. Hence we can restrict ourselves to 
optimize over density matrices which are invariant under the 
action of the symmetry group. 



8 



However, this does not imply that we can restrict the opti- 
mization over pure states (O to pure states with this symmetry, 
since the optimum of the convex roof for a symmetric states p 
is not necessarily taken for symmetric pure states. To see this, 
assume that the problem © gives the correct minimum of the 
problem (f3) for a state IV'min)- Then the symmetric mixed 
state proin := j S/=i QjdV'min}) also satisfies the con- 
straints tr PminWk = Wk and fulfills E(p min ) = E(\ip min )). 
In general, for every IV'min) that minimizes the problem Q 
we find E(p) < E(\ip m i n )) for (without loss of generality) Q 
symmetric states p. Equality holds //fall states in the optimal 
decomposition of p minimize © separately. If we now as- 
sume that a Qj £ Q exists such that Qj l^min) 7^ IV'min), 
i.e. that |^> m in) is not Q symmetric, then the CJ-orbit of all 
minimal pure states constitutes a set of optimal decomposi- 
tion vectors of a flat roof R3l . 

The main problem we encounter in trying to extend the 
above proof for symmetric mixed states to symmetric pure 
states is that the coherent symmetrization J2j ^/QjQj IV'min) 
will typically not satisfy the constraint, unless the witnesses 
have no off-diagonal matrix elements for the Qj \ip m i n ), i.e. 
unless (-0min| Q^WkQi IV'min) oc (5,; ■ If this condition would 



be satisfied, then 



OC Ojl. 

since all decomposition states of p n 



superpositions of the Qj IVViin) - two cases would occur: 



fqjQj iV'min)) = E(t/j min ) for all qj-, V Qj 



a) E(Zi 
1. 



b) A state ipj exists in the range of Q iV'min) such that 

E$) < E(lf) m in). 

From a) we conclude that also the Q symmetric state 
iVsymni) := -jj J2'Li Qj IV'min) is a solution of the 
problem $4^ by virtue of the orthogonality condition 
min) oc Sji . The alternative b) is a contra- 
diction to the initial assumption that |V'min) solves (01). This 
would complete the proof that we can restrict ourselves to op- 
timize over G-invariant pure states. However, the orthogonal- 
ity condition we used in the proof seems a rather stringent at 
first sight; furthermore, in order to test it one would need to 
find the IVWn) first, which amounts to solving the problem 
without imposing the symmetry argument. We leave this dis- 
cussion for future research. Curiously, in the cases considered 
in this work, the minizing pure states turned out to be sym- 
metric. 



applied the latter method to the convex-roof extended three- 
tangle of mixed three qubit states. As a first general result we 
show that a solution to the problem is obtained, when the con- 
straints from the experimental knowledge about the quantum 
state are imposed on pure states instead of mixed states. The 
function convex hull of the thus constructed curve gives the 
solution to the full problem. This elegant simplification can 
be seen as a corrolary to the results in Refs. |@, Si, that has 
not been noticed before. Furthermore, this result highlights a 
close relation to lower bounds obtained from so-called char- 
acteristic curves l22ll . 

We applied both the full and pure state approach to a sim- 
plified situation, where the state is assumed to be of the 
form ( [Tol l, in order to explicitly demonstrate the Legendre 
transformation at work, using three different witnesses. This 
situation corresponds to an experiment with dominating sys- 
tematic errors: in our example, the only ocurring error con- 
sists in the admixture of a W state. This analysis is useful, 
because it also clearly points out the need for a convexifica- 
tion after the related pure-state problem (O has been solved. 
Interestingly, the thus obtained bound coincides with the ex- 
act value of the three-tangle as obtained in Ref. 112811 . Note that 
this restricted approach can be generally applied if a bound on 
the entanglement of a family of states is desired. In addition, 
we have analyzed off-diagonal elements in the witness for the 
same setting, in order to mimick the effect of not properly 
chosen witness operators. The main effect is to reduce the es- 
timated lower bound as compared to a perfect witness, even 
for a pure state at hand. 

We then considered the experimentally relevant case where 
the fidelity of the produced state with respect to the GHZ state 
has been measured, and no prior information about the state 
is assumed to be available. The result is plotted in Fig. [6] and 
has been used to estimate the three-tangle produced in recent 
experiments, summarized in Table I. We only considered ex- 
periments where the fidelity (GHZ| p cxp |GHZ) exceeds 3/4, 
since a lower fidelity is compatible with a state with no three- 
tangle. As to be expected, the witness bound 7 = 5/7 for 
having three-tangle in a GHZ-identity mixture, is reproduced. 

It is a curious observation that the lower bound seems to be 
obtained using permutation symmetric pure states only. Al- 
though the effect of common symmetries of witness and en- 
tanglement measure - as we have here - is directly reflected in 
mixed state solutions of the optimization problem, the same 
conclusion does not seem to be straight forward for pure state 
solutions. 



V. CONCLUSIONS 



We have analyzed the method proposed in Refs. 10, [H for 
obtaining a lower bound for entanglement measures from ex- 
pectation values of witness operators. To this end, we have 
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